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Abstract 

Taking the two-dimensional </> 4 theory as an example, we numerically solve 
the deterministic equations of motion with random initial states. Short-time 
behavior of the solutions is systematically investigated. Assuming that the 
solutions generate a microcanonical ensemble of the system, we demonstrate 
that the second order phase transition point can be determined already from 
the short-time dynamic behavior. Initial increase of the magnetization and 
critical slowing down are observed. The dynamic critical exponent z, the new 
exponent 6 and the static exponents (3 and v are estimated. Interestingly, 
the deterministic dynamics with random initial states is in a same dynamic 
universality class of Monte Carlo dynamics. 

PACS: 05.20.-y, 02.60.Cb, 64.60.Ht, ll.10.-z 



1 



Typeset using REVTJ5X 



It is believed that statistical mechanics is originated from the fundamental equations of 
motion for many body systems or field theories, even though up to now there exists not 
a general proof. For classical systems, equations of motion are deterministic. Statistical 
ensembles are expected to be effective description of the systems. For quantum systems 
the situation is similar. In any case, practically it remains open whether solutions of the 
fundamental equations of motion really yield the same results as statistical mechanics, e.g., 
see Refs. |TJ-§. 

With recent development of computers, to solve numerically equations of motion gradu- 
ally becomes possible. This has attracted much attention of scientists in different areas. For 
example, recently the O(N) vector model and XY model have been numerically solved |6] -Q. 
Noting that energy is conserved during time evolution, solutions of equations of motion ac- 
tually generate a microcanonical ensemble. Making the time average of the observables, 
introducing standard techniques developed in statistical mechanics for the canonical ensem- 
ble, phase transition points and critical exponents have been estimated. The results are in 
agreement with those obtained from a canonical ensemble in statistical mechanics. 

In principle, the deterministic equations of motion should describe not only equilibrium 
properties but also dynamic properties of the systems. In statistical mechanics, dynamics is 
approximately given by some effective stochastic equations of motion, typically at mesoscopic 
level, e.g. Langevin equations or Monte Carlo algorithms. For stochastic dynamics, critical 
slowing down and dynamic scaling around the critical point are characteristic properties. It 
has long been challenging whether stochastic dynamics describes correctly the real physical 
world. It is important and interesting to study dynamic properties of the fundamental deter- 
ministic equations of motion and whether deterministic dynamics and stochastic dynamics 
are in a same universality class. 

On the other hand, in recent years much progress has been achieved in stochastic dy- 
namics. For a long time it is known that around the critical point there exists universal 
scaling behavior in the long-time regime of the dynamic evolution. This is more or less 
due to the divergent correlation time induced by the divergent spatial correlation length. 
Short-time behavior is believed to depend on microscopic detail. However, in recent years it 
is discovered that universal scaling behavior emerges already in the macroscopic short-time 
regime, after a time scale t mic which is sufficiently large in the microscopic sense |p|-[T5|] (for a 
review, see Ref. [fT6|1 ). Important is that one should take care of the macroscopic initial con- 
ditions carefully. For example, dynamic magnetic systems with an initial state of very high 
temperature and small magnetization, the short-time dynamic scaling for the kth moment 
of the magnetization is written as 

M (k \t,T,L,m ) = rW"M (fc) (r 2 i,& 1/, 'r i r 1 L,Pmo). (1) 

Here t = (T — T c )/T c is the reduced temperature, (3, v are the well known static critical 
exponents and z is the dynamic exponent, while the new independent exponent x is the 
scaling dimension of the initial magnetization mo- The parameter b is an arbitrary rescaling 
factor and L is the lattice size. For the dynamics of model A, a prominent property of the 
short-time dynamic scaling is that the magnetization undergoes a critical initial increase at 
early time iH©!! 

M(t)~m t 9 (2) 
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where the exponent 6 is related to xq by 6 = (x — (3/v)/z. 

More interesting and important is that the static exponents /3, v and the dynamic ex- 
ponent z in the short-time scaling form (ED) take the same values as they are defined in 
equilibrium or the long-time regime. This fact makes it possible to determine all these 
exponents already in the short-time regime of the dynamic evolution ||16||. Since now the 



measurements are carried out at early time, the method is free of critical slowing down. The 
short-time dynamic scaling is not only conceptually interesting but also practically impor- 
tant. To clarify whether the short-time dynamic scaling is a fundamental phenomenon of 
the real physical world becomes urgent. 

The purpose of this letter is to study whether there exists short-time universal scaling 
behavior in the dynamics described by the deterministic equations of motion and meanwhile 
determine the critical point and all the static and dynamic exponents. 

The model we choose is the two-dimensional 4 theory. The reason is that the statics 
of this model is known in a same universality class of the two-dimensional Ising model and 
there exist already some numerical results |J. The Hamiltonian of the model on a square 
lattice is 



) 2 - \rn^ + 



(3) 



With TTi 



and it leads to the equations of motion 



E(^+M + ~ 2( t>i) + ™V< - ■ 



(4) 



In the dynamic evolution governed by equations (|j), energy is conserved. The solutions 
are assumed to generate a microcanonical ensemble. The temperature could not be intro- 
duced externally as in a canonical ensemble, but could only be defined internally as the 
averaged kinetic energy. In our short-time dynamic approach, the total energy is actually a 
even more convenient controlling parameter of the system, since it is conserved and can be 
input from the initial state. Therefore, from now r in the scaling form ([!]) will be understood 
as a reduced energy density (e — e c )/e c . Here e c is the critical energy density corresponding 
to a second order phase transition. 

The order parameter of the 4 theory is the magnetization. The time-dependent mag- 
netization M = M^'(t) and its second moment are defined as 



M«(t) 



L 2k 



(k) 



1,2. 



(5) 



The average for observables is only over initial configurations. In the short-time regime 
of the dynamic evolution, spatial correlation length is small, one can easily realize that 
M^ 2) ~ 1/L d . From the scaling form ([!]), we obtain at the critical point for m = ||T~0| , [T8|JT~6 



M {2)^^ t {d-2(}/v)/z_ 



(6) 



Another interesting observable is the auto-correlation function 
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A(t) = ±J2 < u^Ht) > ■ (7) 

i 

From the calculations based on Langevin equations, at the critical point the auto-correlation 



A(t) with m = presents a power law behavior 

A(t) ~ r d/z+e . (8) 

The power law behavior indicates critical slowing down. The interesting point is that the 
new exponent 9 describing the initial increase of the magnetization also appears here. 

Our strategy is that from initial increase of the magnetization, Eq. (Q), we measure the 
exponent 9. Taking 9 as an input, from the power law decay of the auto-correlation, Eq. ([|), 
we extract the dynamic exponent z. Then from power law behavior of the second moment, 
Eq. (Q, we obtain the static exponent (3 jv. Finally, the critical point and the remaining 
exponent v are determined as follows. Assuming L is sufficiently large and t x °^ z is small 
enough, one can deduce from scaling form ([I]) 

M(t) = m t e F(t 1/,/z r). (9) 

When t = 0, Eq. (fj) is recovered. When r ^ 0, the power law behavior will be modified. 
Therefore, searching for a energy density which gives the best power law behavior for the 



magnetization, one determines the critical point |L6 |. The exponent 1/uz can be extracted 
from the logarithmic derivative [TjJ 

d T In M(t,T,m )\ T=0 ~t 1/uz . (10) 

For stochastic dynamics in a canonical ensemble, to obtain the scaling form ([I]) and (||) 
the initial state is taken to be a very high temperature state. Therefore, for deterministic 
dynamics we consider a random initial state with zero or small magnetization. For simplicity, 
we also set initial kinetic energy to be zero, i.e. <pi(0) = 0. Practically we proceed as follows. 
We first fix the magnitude of the initial field to be a constant c, |</>j(0)| = c, then randomly 
give the sign to 0^(0) with the restriction of a fixed magnetization in unit of c, and finally 
the constant c is determined by the given energy. Following Ref. ||, we take parameters 
m 2 = 2. and A = 0.6. 

To solve these equations of motion numerically, we simply discretize <pi by (4>i(t + At) + 
4>i{t — At) — 2<pi{t))/ (At) 2 . After an initial configuration is prepared, we update the equations 
of motion until t = 500. Then we repeat the procedure with another initial configuration. 
In our calculations, we use a fairly large lattice L = 256 and samples of initial configurations 
for average range from 4 500 to 10 000 depending on initial magnetization m®. Errors are 
simply estimated by dividing total samples into two or three groups. Extra calculations with 
L = 128 show that the finite size effect for L = 256 is already negligible small. 

Compared with numerical solutions in the long-time regime, the finite At effect in the 
short-time dynamic approach is less severe, since our updating time is limited. Our results 
are presented with At = 0.05. We have also performed some simulation with At = 0.01 
and confirmed that the finite At effect has been rather small. In any case, however, since 
we have this small At, the calculations are much more lengthy than standard Monte Carlo 
simulations. Therefore, the purpose of this paper is to explore new physics rather than to 
pursue high accuracy of measurements. 
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In order to determine the critical point, we perform simulations with non-zero initial 
magnetization m for a couple of energy densities in the critical regime. In Fig. [l], time 
evolution of the magnetization with m = 0.015 for parameters m 2 = 2. and A = 0.6 has 
been plotted with solid lines for three energy densities e = 20.7, 21.1 and 21.5 in log-log 
scale. To show the universal behavior clearly, we have cut the data for t < t m i C ~ 50. Even 
looking by eyes, one could realize that the magnetization of e = 21.1 has rather good power 
law behavior. Actually, careful analysis of the data between t = 50 and 500 leads to the 
critical energy density e c = 21.11(3). This agrees well with e c = 21.1 given by the Binder 
cumulant in equilibrium in Ref. |J. 

At e c , one measures the exponent 9 = 0.146(3). Accurately speaking, however, the 
exponent 9 is defined in the limit mo — > 0. In general, for finite itlq the exponent 9 may 



show some mo-dependence [ ISf . Therefore, we have performed another simulation with 
m = 0.009 at e c . The magnetization is also displayed in Fig. [3] with a dashed line. The 
corresponding exponent is 9 = 0.158(2). If we linearly extrapolate the results to m = 0, we 
obtain the final value 9 = 0.176(7). 

With the critical energy in hand, we set mo = and proceed to measure the auto- 
correlation A{t) and second moment M^ 2 \t) . In Fig. [| A{t) and M^if) are displayed 
in log-log scale with a solid and dashed line respectively. In order to show some typical 
microscopic behavior within £ TOic , here we present the data from a relatively early time 
t — 10. From the figure we see clearly that at very early time, the curves do not have 
power law behavior, typically oscillating a bit. However, after t > t m i C pa 50, both curves 
nicely stabilize to power law behavior. From the data for t > 50, we measure the exponents 
d/ z - 9 = 0.755(5) and (d - 2/3/v)/z = 0.819(12). 

Finally, to complete the estimate of the exponents, we calculate approximately the deriva- 
tive of the magnetization with respect to the energy density with the data for Fig. [TJ In 
Fig. D, the logarithmic derivative of the magnetization is plotted in log-log scale at the crit- 
ical point. As usual, the fluctuation here is larger The power law behavior is also not 
as clean as that of the previous observables. The slope of the curve tends to be smaller 
as the time t evolves. To improve this situation, probably we need to perform simulations 
with energy densities closer to e c and with high statistics and/or to consider corrections 
to scaling. This requests very much computer times. In any case, from the slope between 
t = 200 and 500, we estimate 1/uz = 0.492(26). 

In Table [IJ we summarize all the values of the critical exponents discussed above. From 
these values, we can estimate the critical exponents z, 2(3 /v and v. The results are also 
listed in Table |I[ For comparison, the exact values of 2f3/v and v for the Ising model and 
available results for other exponents measured from a similar dynamic process with standard 
Monte Carlo dynamics are also given in Table |. Remarkably, not only the static exponents, 
but also the dynamic exponents of the 4 theory are very close to those of the Ising model 
with standard Monte Carlo dynamics. The exponent 9 shows some percents of difference. 
However, taking into account that the measurements of 9 are not so easy, especially the 
extrapolation to mo = limit may induce some systematic errors, we would consider that 
the exponent 9 for both model are the same. Therefore, we conclude that the 4 theory 
described by the deterministic equations of motion are in a same static as well as dynamic 
universality class of the Ising model with standard Monte Carlo dynamics. 

Why is the deterministic dynamics in a same universality class of the stochastic dynam- 



5 



ics? The origin may be traced back to the random initial state and the kinetic energy term 
in the Hamiltonian. The random initial configuration of {0j} induced stochastically evolv- 
ing kinetic energy. The stochastic kinetic energy serves as a kind of noises or heat bath to 
the potential energy. This is similar to the stochastic dynamics with a canonical ensemble, 
where the Hamiltonian is simply the potential energy term here. 

In conclusions, we have numerically solved the deterministic equations of motion for the 
two-dimensional 4 theory. Short-time universal behavior is found. Based on the short-time 
dynamic scaling form, the critical point and all the static and dynamic critical exponents are 
determined. The values of the static exponents agree with those of the Ising model. More 
interestingly, both the dynamic exponents z and 6 also coincide with those of the kinetic 
Ising model induced by the standard Monte Carlo algorithms (heat-bath and Metropolis 
et al). Especially, the dynamic exponent z = 2.15(2) is very far from z — 1, which is the 
expectation from the naive deterministic viewpoint. These results show that, on the one 
hand, the deterministic equations of motion indeed can describe statistical properties of the 
systems both in statics and dynamics, on the other hand, the Monte Carlo dynamics can be 
a good effective description of the real physical dynamics, at least in some cases. 

Since the measurements are carried out in the short-time regime of the dynamic evolu- 
tion, our dynamic approach does not suffer from critical slowing down. Furthermore, the 
errors induced by a finite At in our measurements are also limited. It is challenging to 
derive analytically the short-time dynamic scaling from deterministic equations of motion. 
Extension of the present work to quantum systems would be very interesting. 

Acknowledgements: Work supported in part by the Deutsche Forschungsgemeinschaft; 
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FIGURES 
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FIG. 1. Time evolution of the magnetization in log-log scale. Solid lines are for mo = 0.015 
with energy densities e = 20.7, 21.1 and 21.5 (from above), while the dashed line is for mo = 0.009 



with e c = 21.11. 
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FIG. 2. The auto-correlation A(t) (solid line) and second moment M^(t) (dashed line) with 
mo = at the critical point in log-log scale. 
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FIG. 3. Logarithmic derivative of the magnetization at the critical point in log-log scale. 
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TABLES 





e 


d/z-e 


(d-2/3/u)/z 


l/vz 


z 


2/3/u 


V 




0.176(7) 


0.755(5) 


0.819(12) 


0.492(26) 


2.148(20) 


0.24(3) 


0.95(5) 


Ising 


0.191(1) 


0.737(1) 


0.817(7) 




2.155(3) 


1/4 


1 



TABLE I. The critical exponents measured for the cj) theory in comparison with those of the 
Ising model. The values of and v for the Ising model are exact, while others are taken from 
Table 2 in Ref. ||. 
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